Pulsations in Short GRBs from Black Hole-Neutron Star Mergers 



Nicholas Stone, Abraham Loeb, and Edo Berger 
Astronomy Department, Harvard University, 60 Garden St., Cambridge, MA 02138, USA 

(Dated: September 20, 2012) 

The precise origin of short gamma ray bursts (SGRBs) remains an important open question in 
relativistic astrophysics. Increasingly, observational evidence suggests the merger of a binary com- 
pact object system as the source for most SGRBs, but it is currently unclear how to distinguish 
observationally between a binary neutron star progenitor and a black hole-neutron star progeni- 
tor. We suggest the quasi-periodic signal of jet precession as an observational signature of SGRBs 
originating in mixed binary systems, and quantify both the fraction of mixed binaries capable of pro- 
ducing SGRBs, and the distributions of precession amplitudes and periods. The difficulty inherent 
in disrupting a neutron star outside the horizon of a stellar-mass black hole biases the jet precession 
signal towards low amplitude and high frequency. Precession periods of ~ 0.01 — 0.1 s and disk-BH 
spin misalignments ~ 10° are generally expected, although sufficiently high viscosity may prevent 
the accumulation of multiple precession periods during the SGRB. The precessing jet will naturally 
cover a larger solid angle in the sky than would standard SGRB jets, enhancing observability for 
both prompt emission and optical afterglows. 



I. INTRODUCTION 

The origin of the short-hard gamma-ray bursts 
(SGRBs; durations T go < 2 s; pQ) was largely a mat- 
ter of speculation until the recent discovery of their af- 
terglows and host galaxies (e.g., [2H5]). These observa- 
tions have demonstrated that SGRBs are cosmological 
in origin (z > 0.1; [5]); have a beaming-corrected energy 
scale of ~ 10 49 — 10 50 erg [ZHS]; lack associated super- 
novae (e.g., [SI I10J ) ; occur in a mix of star- forming and 
elliptical galaxies |llj ; have a broad spatial distribution 
around their hosts [12] , with some events offset by tens of 
kpc [13] ; and have low-density parsec-scale environments 
[51 [T2] . The confluence of these characteristics provides 
support to the popular model of compact object (CO) 
mergers [14) . 

In this context, the key open question which motivates 
our paper is: if SGRBs originate in CO mergers, what 
types of compact objects are merging? Specifically, a 
neutron star (NS) can theoretically be tidally disrupted 
by, and produce an accretion disk around, either a more 
compact NS, or a sufficiently small stellar mass black 
hole (BH). Unfortunately, it is not clear how to distin- 
guish between mergers of neutron star binaries (NS-NS) 
and mixed binaries (BH-NS). The advent of gravitational 
wave (GW) astronomy will facilitate this task, since the 
waveforms accompanying NS-NS mergers may be suffi- 
ciently distinct from those in BH-NS mergers, but the 
Advanced LIGO era is still half a decade away. It is also 
possible that GW signals seen by Advanced LIGO will 
lack accompanying electromagnetic counterparts if the 
SGRB beaming angle is too small, or the intrinsic event 
rate too low. In this paper we suggest a clear observa- 
tional tool for distinguishing between NS-NS and BH-NS 
progenitors of SGRBs using electromagnetic data only: 
the precession of the disks and jets associated with these 
events. 

Jet precession has previously been discussed as a phe- 
nomenon relevant for GRBs, originally with regard to 



now-disfavored GRB models [TBI HI] but later in the 
context of CO mergers. Early works in the CO merger 
paradigm considered disks fed by stable mass transfer 
from a NS onto a BH so that the precession was forced 
by tidal torques [IB], but subsequent models considered 
the more realistic neutrino-dominated accretion flows 
(NDAFs) formed after tidal disruption of a NS by a stel- 
lar mass BH [Jj5]. This last model is similar to the one 
presented in this paper, in that it considers a thick disk 
precessing as a solid body rotator due to general relativis- 
tic Lense-Thirring torques. It has been applied to pre- 
dict different observational signatures, such as the light 
curves of precessing jets [JJJ][2U] or even LIGO-band grav- 
itational wave (GW) signals emitted by forced precession 
of a large amount of disk mass [5TJ [23] . Others have con- 
sidered an inclined disk whose inner regions include a 
Bardeen-Petterson warp; in their model the inner region 
of a neutrino-dominated accretion flow (NDAF) precesses 
along with the BH spin about the total angular momen- 
tum vector [25] . 

This paper follows Ref. [lj5] in focusing on thick 
disks precessing as solid body rotators, which are well- 
motivated for the supercritical accretion flows and mis- 
aligned angular momentum vectors characteristic of com- 
pact object mergers (§2). Our work differs from past ef- 
forts, however, in our consideration of the viscous spread- 
ing of the disk, as well as our adoption of simplifying as- 
sumptions tailored to match results from numerical rel- 
ativity (NR) simulations of mixed binary mergers (§3). 
We quantify for the first time the distributions of pre- 
cession periods and angles given physically motivated as- 
sumptions about progenitor spins and masses (§4), us- 
ing simple analytic formulae when appropriate and more 
complex empirical fits to NR simulations when necessary. 
In the process we estimate the fraction of BH-NS mergers 
which can actually produce accretion disks. We conclude 
by considering the observable consequences of jet preces- 
sion in the context of SGRBs (§5). Unlike in previous 
work, we discuss both the case where the jet is tied to 



the BH spin vector, and the case where it aligns with the 
disk angular momentum vector. 



II. DISK PRECESSION 

Non-axisymmetric torques will, initially, induce small 
warps in accretion disks due to differential precession be- 
tween adjacent mass annuli. The evolution of these warps 
depends on how they are able to propagate through the 
disk. When the timescale for precession is shorter than 
the sound crossing time, these warps will propagate in a 
diffusive way, allowing differential precession to produce 
substantial shear viscosity and dissipating large amounts 
of orbital energy. In the context of tilted accretion disks 
around spinning black holes experiencing Lense-Thirring 
torques, diffusive propagation of warps will align the in- 
ner regions of the disk with the black hole equatorial 
plane; this is known as the Bardeen-Petterson effect [2~11 - 

1251 . 

In the opposite regime, a thick disk with a short sound- 
crossing timescale will propagate warps in a wavelike 
manner, redistributing torques throughout the disk and 
inducing near-rigid body rotation [27l|28]. In particular, 
rigid body rotation is possible if H/r > a, where H/r is 
the disk height and a the dimensionless viscosity param- 
eter at a radius r. This is the regime most relevant for 
compact object mergers, and is therefore what we will 
consider for the remainder of this paper. Approximately 
rigid body precession has been seen in hydrodynamical 
simulations of protoplanetary disks being torqued by a 
binary companion [29] , in GRMHD simulations of tilted 
accretion disks around spinning black holes [301 |3I] , and, 
notably, in NR simulations of BH-NS mergers [32] . 

In the Newtonian limit, a solid body rotator will pre- 
cess with a period T proc = 2-7T sin ipd(J / T ), where ipd is the 
misalignment angle between the accretion disk and the 
BH equatorial plane, J is the total angular momentum 
of the disk, and r is the Lense-Thirring torque integrated 
over the entire disk. Specifically, if the disk possesses a 
surface density profile that is nonzero from inner ra- 
dius Ri to outer radius R , and the disk elements possess 
orbital frequency f2(r), then 

f Ro 

J = 2ir S(r)fi(r)r 3 dr, (I) 

and 

T = 47r G2M BH«BH s . n ^ f R ° JL E (r)fi(r)r a dr, (2) 

c A J R . r 6 

where the BH's mass and dimensionless spin are Mbh 
and obh, respectively. Throughout this paper G is the 
gravitational constant and c is the speed of l ight. If we 
use the Keplerian orbital frequency £lk = \J GA/bh/?' 3 , 
then for a density profile of the form S(r) = £o(r/ro)^, 
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the precession timescale is 

^ 7rr s (l + 2Q r° /2 -S 1/2+C (l - (n/r.) 5 ^) 
proc c(5-2<) a B n(l-(ri/r )i/ 2 +t) ' 

(3) 

Here we have normalized r Q — R /r g and ri = R[/r s by 
the gravitational radius r g = GA/bh/c 2 - Other effects 
which influence Eq. ^ include nutation and relativistic 
corrections to the orbital frequency, which we neglect in 
this analysis (but see [T51 119| ). 

We stress that T prec is an instantaneous precession pe- 
riod, and that for the non-steady state disks expected in 
CO mergers important quantities such as £ and r Q will be 
time-dependent (although the cancellation of E in J/t 
means that the secular decrease in disk mass will not af- 
fect T prcc ). The time dependence of these variables will 
cause any signal to be quasi-periodic rather than peri- 
odic. Because the dominant feature of the disk's evolu- 
tion will be viscous outward spreading [331 134] , we expect 
Tprcc to increase with time. 

Although a power law definition of £ is appealing for 
its simplicity, both analytical models [33 and NR simula- 
tions [32] indicate that the true structure of these disks is 
more complex. To better account for realistic disk struc- 
ture, and also to quantify the time evolution of the disk 
as it spreads outward, we adopt the SGRB disk model 
of Ref. [31], which derives exact £ solutions for a vis- 
cously spreading ring of matter, and then couples these 
solutions to more detailed models of disk energetics and 
composition. As matter from the disk accretes onto the 
BH, the intially advective, neutrino-dominated accretion 
flow will become optically thin to neutrinos and geometri- 
cally thinner. At later times (t > 0.1 sec) the disk will be- 
come a geometrically thick, radiatively-inefficient accre- 
tion flow (RIAF). In principle, the intermediate neutrino- 
cooled period could prevent later disk precession by align- 
ing the disk into the BH midplane through the creation 
of a Bardeen-Petterson warp. In practice, it seems that 
even the neutrino-cooled phase of accretion still possesses 
H/r > a, and is therefore unlikely to align (Ref. [33] . 
Fig. 2). 

We are primarily interested in the RIAF stage, both 
because it has the longest duration, and because for low 
disk masses (Md < O.IMns, which is the case for most 
BH-NS mergers - see §111), it is the only phase of accre- 
tion. For this stage of disk evolution, we can write the 
surface density as 

(2x 1 ~ n / 2 \ 

X h/\i-2n\ I " I ■ 

Here Md is the initial disk mass, R^is is the initial radius 
of the spreading mass ring (i.e. the radius where the NS 
is disrupted), I m is a modified Bessel function of order 
to, x = r/R dis , t = t(12iv (l — n/2) 2 / R\ is ), and we have 
assumed viscosity of the form v = vqx ti . We calibrate Vq 
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with the initial relation t v isc,o = R\i s l v ano - ^ ne equation 



*visc,0 ~ 6 x 10" 



0.1 



3M a 



10 5 



(5) 



0.3R, 



(lis 



where a is the dimensionless Shakura-Sunyaev viscosity 
coefficient and Ho is the characteristic disk height. The 
value of a is set by the magneto-rotational instability 
(MM), and has been estimated to span a wide range of 
values, from ~ 0.01 in local, shearing box simulations 
[33] to ~ 1 in global GRMHD simulations [35] • However, 
large a values seen in global simulations are confined to 
small radii, and in these simulations a ~ 0.1 at r > 10r g . 
The importance of a for our results lies primarily in how 
viscosity controls the outward spreading of the disk, so 
we follow Ref. [37] and consider large-radii a values of 
0.01, 0.03, and 0.1. Hq will vary both in radius and in 
time; generally, Hq grows as one moves further out in the 
disk [37] , and also as the outer edge of the disk viscously 
spreads, putting a larger fraction of the disk into a purely 
advective regime with large height [34) . Our results are 
fairly sensitive to both a and Hq, but because the for- 
mer spans a wider range we vary a and fix Hq = 0.3-R. 
We arrive at this value by considering the size of the 



disk at a time ti/2, a characteristic, "halfway," preces- 
sion timescale. Specifically, t 1 / 2 — ((^visc o + l)/^) -3 is 
the time at which half the SGRB's precession cycles will 
have occurred if it lasts from t — i v isc,o to t = 1 s and 
viscous spreading of the disk causes T prec oc t 4 / 3 . Typi- 



cally, ii/2 ~ 100 msec, which corresponds to a disk outer 
edge at r ~ 50; at this distance, more detailed modeling 
of disk structure [37] indicates H/R sa 0.3. 

Using Eq. Q, we plot the time evolution of T prcc in 
Figure [T] and find that it increases in rough agreement 
with analytic expectations: at late times, Eq. Q ap- 
proaches a power law with Q — 1/2, and the outer edge 
of the disk r oc i 2 / 3 . Assuming that r; remains fixed (and 
ignoring lower-order contributions from r Q ), Eq. ^ then 
implies T proc oc t 4 / 3 . In Fig. Ill we also plot N ( 

rcle 



* cycles ■ 



the 



total number of precession cycles undergone during the 
GRB. For a > 0.1, a viscously spreading SGRB disk will 
generally experience N cyc \ es < 1. 

In the above discussion we have assumed that the an- 
gular momentum lost by inspiralling disk matter is re- 
distributed outward by internal viscous torques. An al- 
ternate possibility is angular momentum loss through a 
magnetized disk wind [33], in which case r oc i 2 / 5 and 



T, 



oc t 4 / 5 . If this is the case, N cyc i cs > 10 for all real- 
istic a values. For the remainder of this paper, however, 
we conservatively calculate fiducial precession timescales 
using Eq. Q and n = 1/2. Our results are generally in- 
sensitive to i?dis and independent of (although we will 
later use as a criterion for whether or not an SGRB 
can form - see §111). On the other hand, our results are 
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FIG. 1. Panel (a) shows the time evolution of T proc assuming a 
viscously spreading disk structure given by Eq. Q. The black 
dotted curves represent a = 0.1, the dashed magenta curves 
a = 0.03, and the solid cyan curves a — 0.01. Thick curves are 
for nearly equatorial disruptions with obh = 0.9, while thin 
curves are for <xbh = 0.9 and initial spin-orbit misalignment of 
70°, or equivalent!/ a nearly aligned disruption with a ~ 0.5. 
The dash-dotted red line is oc t A ^ , the rough time evolution 
of Tprec- Panels (b) and (c ) show AT cyc i es , the accumulated 
number of cycles for 0.1 s < f < 1 s and 0.01 s < t < 1 s, 
respectively. 



fairly sensitive to a, with large a values increasing the 
precession period and decreasing the number of preces- 
sion cycles that can fit in the duration of an SGRB. In 
all cases we self-consistently calculate the inner disk edge 
n using the formalism in Perez-Giz 2012 |38. for finding 
the innermost stable spherical orbit (ISSO), the tilted 
analogue to the innermost stable circular orbit (ISCO). 
Details of the ISSO calculation are in Appendix A. 

The precession of the SGRB disk in isolation is unlikely 
to be observable, and is mainly interesting as a source of 
jet precession. The observational signatures of jet preces- 
sion will hinge on two uncertain astrophysical questions: 
the opening angles of SGRB jets, and the alignment di- 
rection of a jet in a tilted accretion flow. The first of these 
questions has recently become amenable to observational 
constraint; observations of jet breaks in SGRBs suggest 
opening angles of ~ 10° [5]. Observational evidence for 
the second question is limited, and ambiguous. Observa- 
tions of a relativistic outflow following tidal disruption of 
a star by a supermassive BH (Swift J164449.3+573451) 
suggested that in that case, the jet aligned with the BH 
spin axis |39j ; on the other hand, observations of the mi- 
croquasar LSI+61303 have been interpreted as evidence 
of a precessing jet, aligned with the angular momentum 
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axis of a processing disk [3tJ]. There may not be a uni- 
versal answer to this question, as different hypothetical 
jet launching mechanisms might each tie the jet axis to 
a different preferred direction. 

However, for the two leading jet launching mechanism 
candidates in SGRBs - vv pair annihilation [UH2], ano - 
the Blandford-Znajek (BZ) mechanism [32 EEH - there 
are theoretical reasons to believe that the jet will align 
with the disk angular momentum vector. The vv anni- 
hilation scenario is independent of BH spin and depends 
only on disk properties. Alignment of a BZ-powered jet 
is more ambiguous, but recent works that have consid- 
ered jet precession in SGRBs assumed that a jet powered 
by the BZ mechanism will align with the disk angular 
momentum vector [Till [211] because the magnetic field is 
anchored in the disk. This assumption has been further 
supported by NR simulations of force-free electromag- 
netic fields around spinning BHs |45| . which found that 
the direction of Poynting flux from the BZ mechanism 
is governed by larger-scale magnetic fields and not the 
BH spin vector, although we note that GRMHD simu- 
lations of tilted accretion flows with matter have so far 
been unable to resolve jets [46] . 



III. PROGENITOR BINARIES 

Because large amplitude precession requires large am- 
plitude misalignment of the post-merger BH and its ac- 
cretion disk, we must consider which CO mergers can 
actually produce misaligned disks. Despite their larger 
pre-merger spin-orbit misalignment, NS-NS mergers are 
unlikely to produce significantly misaligned disks. If one 
or both members of the NS-NS binary were millisecond 
pulsars, disk precession could be feasible: spin angular 
momentum Jmsp ~ §-^ns-Rns^ns ~ 1 x 10 42 kg m 2 /s, 
and the orbital angular momentum at the disruption ra- 
dius I/NS-NS 

10 43 kg m 2 /s. Thus, the misalignment angle ip^ between 
the post-merger accretion disks (which we assume to lie 
in the initial orbital plane) and the BH spin axis will be 
~ 5° in NS-NS mergers involving one star with a ms spin 
period, if we assume initially orthogonal spin and orbital 
angular momentum vectors. Standard population syn- 
thesis channels indicate, however, that the one recycled 
component of NS binaries typically has a minimum spin 
period of w 4 ms [47) . which would imply ip^ ^5 1°, a value 
that is likely too small to carry significant observational 
consequences. 

In BH-NS mergers, however, the BH may possess a 
larger natal reservoir of spin angular momentum, allow- 
ing for greater misalignment between the post-merger 
BH and the disk formed from NS debris (which wc 
have assumed to lie in the initial orbital plane). Na- 
tal spin is the most relevant quantity, although subse- 
quent mass transfer onto the BH may produce modest 
changes in obh [3H]- For a BH-NS system, the rele- 
vant numbers are Jbh = «bhGMbh/c and Lbh-ns ~ 



(1 + q)M BH {R NS ) V(l + g)GM BH /(i?Ns) 3 , where a BH 
is the dimensionless black hole spin and we define the 
mass ratio q = M NS /M BH - For M BH = 5M (1OM ) and 
flBH = 0.9, the disk misalignment angle V>d ^ 20° (30°), 
which is large enough to be observationally interesting. 
Here we have also assumed initially orthogonal spin and 
orbital angular momentum vectors. 

These simple Newtonian estimates motivate an inves- 
tigation of BH-disk misalignment in BH-NS mergers, 
but are insufficient for accurately estimating either post- 
merger BH spin a BH or the misalignment angle be- 
tween a BH and L B h-ns- This is because they do not 
account for the fully dynamical, strong field GR effects 
that accompany compact object mergers, and as we shall 
see underestimate post-merger misalignment. Recently, 
empirical post-Newtonian (PN) formulas were derived to 
calculate these quantities in the case of BH-BH mergers 
[50) . These formulas have a large number of free param- 
eters which were calibrated based on a suite of numerical 
relativity (NR) simulations. Because disruptions of NSs 
by BHs are so marginal (i.e. occur so close to the ISSO), 
it is reasonable to expect these formulas to have some 
utility in making predictions for BH-NS mergers; we show 
in Appendix B that they are in fact quite accurate when 
compared to NR simulations of BH-NS coalescence. For 
this reason, we use these PN formulae to calculate the 
final spin a BH for a BH formed by the merger of a BH- 
NS with mass ratio q, pre-merger BH spin a B n, and a 
pre-merger spin-orbit misalignment angle ip. We approx- 
imate the pre-merger NS spin magnitude as <3ns = 0. Fi- 
nally, we calculate cos ip d = o BB ■ iBH— ns/ (cbh-^bh— ns) • 

Even though the PN fomulae in Appendix B repro- 
duce a BH with reasonable accuracy, our assumption that 
the post-merger disk tilt V>d is described by the angle be- 
tween the initial orbital plane and a BH may overestimate 
tpd- When we compare our calculations to the late-time 
(40 ms) results of Ref. (32], we appear to overestimate 
■0d by a factor ~ 3, although our approach provides a sig- 
nificantly more accurate estimate of early-time disk tilt. 
Similar evolution of disk tilt is not seen in steady-state 
GRMHD disk simulations [30], but is likely at least par- 
tially physical for BH-NS mergers, as the disk adjusts 
to an equilibrium configuration. However, numerical vis- 
cosity may also play a role in reducing the disk tilt in 
misaligned NR simulations (Francois Foucart, personal 
communication). For simplicity, we do not model the 
time evolution of ^d, but this probably results in overes- 
timates of ~ 2, which we note again in §V. 

The final component in our calculation is a criterion 
for SGRB production in a BH-NS merger. The tidal ra- 
dius, defined in the Newtonian limit as r t = -Rns? _1 ^ 3 j 
is appealing for this purpose: only a fraction /grb of all 
BH-NS mergers will produce an accretion disk and jet, 
because if r t < nsso, the NS is swallowed whole. How- 
ever, while the tidal radius is cleanly defined in other 
contexts [51j . in the case of relativistic, comparable mass 
mergers it is not obvious that the Newtonian definition 
is applicable [52) . Furthermore, the gravitational radius 
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r g ~ i?NS- In light of these complications, we adopt a 
fitting formula for the initial remnant disk mass Ma , cal- 
ibrated from NR simulations of aligned BH-NS mergers 



Ma 



M] 



0.415q 



-1/3 



NS 



NS 



C 2 i?NS 7 



0.148^^. (6) 
Rns 



Although this fitting formula was calibrated from NR 
data on aligned mergers, we generalize it to misaligned 
mergers by substituting nsso f° r nsco- This appears 
to reproduce NR simulations of misaligned postmerger 
disks reasonably well (Appendix B). 

The limiting value of M^/M^s required to produce a 
SGRB is highly uncertain, but past theoretical work as- 
suming jets are powered by vv annihilation has suggested 
that SGRBs are viable for M d /M NS > 0.01 [42]; likewise, 
a recent attempt to observationally infer disk masses as- 
suming vv annhilation [53] found 0.01 < M^/M^s < 0.1. 
In this work, we take Md/M-^g > 0.05 as the cutoff for 
SGRB production, but discuss the effects of stricter and 
weaker criteria in SIV and V. 



IV. DISTRIBUTIONS 

We now integrate the above analytic criteria over dis- 
tributions of progenitor masses and spins to find distribu- 
tions of /grb, tpd, Ancles, and T pTCC (t 1/2 ). Because the 
distributions of progenitor quantities are not at present 
precisely constrained by observation or population syn- 
thesis, we consider a wide range of possibilities to bracket 
the available parameter space. For our fiducial case, we 
take the parametric BH mass function from Ref. [55] 
(hereafter the "OPNM mass function"), given by 



Pi 



OPNM 



(M BH ) 



e M c /M scale 
M sca \ c 



-M/M scale 



0. 



M > M c 

M < M r 



where the best-fit values were found to be Af sca i e ~ 
1.7-Mq and M c w 6.2M & [55]. An important qualitative 
feature of the OPNM mass function is the large mass 
gap between NSs and the lowest-mass BH. Motivated by 
recent results [50] which suggested the mass gap may be 
less distinct than in Ref. [55J , we consider as an alternate 
case a Gaussian mass function where the best-fit values 
for mean mass /xbh and dispersion ctbh were found to be 
7.35M and 1.25M Q , respectively [5Tj . 

Recent observations have measured spins for 7 BHs 
and placed upper or lower limits on spins for 3 more |62| - 
[53] . Although observations of more systems are needed, 
the current spin distribution is noticeably bimodal. Be- 
cause of the small number of data points we do not at- 
tempt to fit a parametrized spin function, and instead 
simply take a flat prior on asHj sampling it uniformly in 
the ranges (0, 0.3) and (0.7, 1) for our fiducial, "bimodal" 
case. For non-fiducial cases, we also consider three alter- 
nate spin functions. The "flat" , "slow," and "fast" cases 



uniformly sample asH along the intervals (0,1), (0,0.5), 
and (0.5, 1), respectively. 

Kicks resulting from asymmetric supernova explosions 
are expected to produce spin-orbit misalignment in BH- 
NS binaries. Past research has constrained the allowed 
pre-merger misalignment angle ip as a function of pro- 
genitor masses and separation, and kick velocity distri- 
butions [3S] . More recent population synthesis of BH-NS 
binaries has found a wide spread in pre-merger spin-orbit 
misalignment ip, but with ~ 50% of systems possessing 
ip < 45° |55] • O ur fiducial, "prograde" case samples the 
pre-merger spin-orbit misalignment uniformly in ip be- 
tween 0° and 90°, but we also consider an alternate, 
"isotropic" case where ip is sampled uniformly from 0° 
to 180°; physically this would represent larger supernova 
kicks. 

We sample NS masses from a Gaussian distribution 
peaked at a mean //ns = 1-35M0 with standard devia- 
tion (Tns = 0.13M Q . These values, taken from the double 
NS binaries examined in Ref. [57], are in good qualita- 
tive agreement with other studies of the NS mass func- 
tion [55]. Because most NS equations of state that are 
not in conflict with observations of ss 2M® NSs [55] are 
roughly constant-radius in the relevant mass range, we 
take a fiducial radius of 13.5 km, but as an alternate case 
consider a NS radius of 11 km. 

With these distributions defined, we are now ready to 
populate a large Monte Carlo sample of BH-NS mergers. 
Our precise procedure is as follows, for any desired set of 
distributions: 

1. Generate 2 x 10 5 BH masses Mbh, spin magnitudes 
<2bh, initial misalignment angles ip, and NS masses 
M N s- 



2. Compute the pre-merger nsso from Eq. (A3 1. 

i Mgjj, BH spin 



3. Calculate the post-merger BH mass -fw tn 

q RH , and spin-disk misalignment ipd using Eqs. 
JB2| and p3| . 



4. Flag the disruption as GRB-producing if 
M d /M NS > 0.05. 

5. Compute iV cyc i cs , T pTec (t 1/2 ), and /grb (using the 
post-merger r' lsso ). 

In Table I, we summarize the scenarios considered in 
this work, along with averaged results. Typical preces- 
sion periods are (T prec (ti/2)) ~ 50 ms — 100 ms, with 
mean disk-BH misalignment typically (tpd) ~ 20°. These 
results are generally insensitive to variation of assump- 
tions about the progenitor population. The fraction of 
BH-NS mergers which produce SGRBs (/grb) varies 
more with the properties of progenitor binaries, and in 
particular depends strongly on the value of pre-merger 
BH spin obh- Values of obh ^ 0.5 are generally neces- 
sary to produce a post-merger disk. Scenario A, our fidu- 
cial case, results in GRBs for sa 30% of BH-NS mergers. 
Results for scenario A are quite similar to those in sce- 
narios B (Gaussian BH mass function) and D (flat prior 



6 



Scenario 


BH Masses 


BH Spins 


i> 


M d /M NS 


Rns 


/grb 


{-^i cycles) 


(7prec(£l/2)) 


(V>d> 


A 


OPNM 


Bimodal 


Prograde 


0.05 


13.5 km 


0.303 


6.3 (14.7,2.2) 


48 (38, 58) ms 


20.9° 


B 


Gaussian 


Bimodal 


Prograde 


0.05 


13.5 km 


0.333 


5.9 (13.7,2.1) 


54 (43, 66) ms 


20.9° 


C 


OPNM 


Bimodal 


Prograde 


0.05 


11 km 


0.131 


7.4 (17.2,2.6) 


34 (27,42) ms 


15.9° 


D 


OPNM 


Flat 


Prograde 


0.05 


13.5 km 


0.229 


5.7 (13.2,2.0) 


57 (45, 69) ms 


19.0° 


E 


OPNM 


Slow 


Prograde 


0.05 


13.5 km 


0.010 


2.6 (5.9,0.9) 


131 (105,157) ms 


9.3° 


F 


OPJNM 


Fast 


Prograde 


0.05 


13.5 km 


0.447 


5.8 (13.4, 2.1) 


55 (44, 67) ms 


19.3° 


G 


OPNM 


Bimodal 


Isotropic 


0.05 


13.5 km 


0.248 


4.8 (11.2,1.7) 


64 (51,78) ms 


32.0° 


H 


OPNM 


Bimodal 


Isotropic 


0.05 


11 km 


0.099 


6.8 (15.8,2.4) 


38 (29,46) ms 


23.2° 


I 


OPNM 


Bimodal 


Prograde 


0.01 


13.5 km 


0.348 


6.1 (14.1,2.2) 


52 (42, 63) ms 


22.6° 


J 


OPNM 


Bimodal 


Prograde 


0.2 


13.5 km 


0.102 


7.7 (18.1,2.7) 


31 (24, 39) ms 


15.3° 



TABLE I. Choices of black hole mass, black hole spin, and spin-orbit misalignment distributions we have surveyed in this 
paper. We also summarize key results here: the fraction of BH-NS mergers that can produce a disk and GRB (/grb), a typical 
precession period in milliseconds (T P rec(ii/2)), an d the mean post-merger misalignment angle between the disk and the BH 



equatorial plane (ipd)- The first number printed in the (T prec } and {N cyc i cs ) columns is the fiducial value for ot 
the numbers in parentheses represent a = 0.01 and a — 0.1 cases. 



0.03, while 




20 40 60 80 20 40 60 80 20 40 60 80 25 50 75 100125150175 

<Ad <Ad <Ad <Ad 



FIG. 2. The probability distribution of post-merger misalignment angles ipd for scenarios A, B, C, D, E, F, G, and H (each 
panel is labeled accordingly). The vertical axis is a probability and the horizontal axis is ipd- The dark blue curves show the 
ipd distribution for all 2 x 10 BH-NS mergers in each scenario, while the smaller, light orange curves show only the subset of 
events where the NS is disrupted outside the ISSO (i.e. the subset where a disk and jet can form). Generally, /grb, the ratio 
of the light orange area to the dark blue area, falls for bottom-heavy spin distributions and smaller NS radii. All distributions 
of ipd cut off fairly sharply above 50° . 



on the BH spin distribution). The probability of disrup- 
tion of the NS falls by a factor w 2 in scenario C (softer 
NS equation of state), and becomes almost negligible in 
scenario E (bottom- heavy spin distribution). The frac- 
tion of BH-NS mergers capable of producing SGRBs is 
maximized in scenario F (fast spins), reduced somewhat 



if we switch to an isotropic ip distribution (scenario G), 
and reduced significantly for a stricter lower cutoff on 
Md/MNS (scenario J). Our fiducial scenario comes close 
to maximizing /grb; only switching over to top-heavy 
BH spin distributions (scenario F) or laxer SGRB crite- 
ria (scenario I) increase its value. More detail on /grb 
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FIG. 3. The probability distribution of post-merger misalign- 
ment angles ipd for scenarios I and J, illustrating the depen- 
dence of /grb and on the strictness of our requirement 
for M(j/Mns. Specifically, in the top panel we have imposed 
the strict requirement that M^/M^s > 0.2 in order to pro- 
duce a SGRB, while in the lower panel we have imposed the 
much laxer requirement of Md/M-^s > 0.01. In all other re- 
spects both of these cases are identical to scenario A. Results 
in the bottom panel are quite similar to scenario A, but in the 
top panel, /grb has been strongly suppressed, particularly at 
higher ipd- The axes and curves are the same as in Fig. [2] 



and distributions of ipd can be seen in Figure [2j 

In Figure |3j we explore the sensitivity of /grb and ipd 
to our uncertain assumption about the minimum normal- 
ized disk mass, Md/M^s, needed to produce a SGRB. In 
our standard scenarios, we take a cutoff value of 0.05, 
but in scenarios I and J we change this cutoff to 0.01 and 
0.2, respectively. While relaxing the cutoff seems to have 
little effect on our overall results, increasing the cutoff 
above Md/M-^s ~ 0.1 very quickly suppresses /grb, and 
biases those SGRBs which are produced toward low-tilt 
accretion disks. 

In general, ipd is strongly cut off above 40° — 50° in all 
scenarios except E, and there is little variation between 
the mean value (^d) in each of our progenitor scenarios. 
Specifically, (ipd) ranges from 9° to 32°. Interestingly, our 
two scenarios with isotropic pre-merger ip only modestly 
extend the range of post-merger ipd values: it seems that 
strongly misaligned (ip > 90°) BH-NS mergers simply 
do not produce massive disks. As mentioned earlier, our 
calculation of ipd likely overestimates its true value by a 
factor « 2, because the tilt angle will decrease as the disk 
settles into equilibrium, so the true range of (ipd) would 
more accurately be ~ 5° to ~ 15°. For rapidly spinning 
BHs, it is possible that our ipd estimate accurately de- 
scribes the first few precession cycles, before disk tilt has 
time to decrease. 

As can be seen in Table I, the key observables (2V cycles ) 
and (T picc (ti/ 2 )) depend much more sensitively on a than 
they do on the parameters of BH-NS binary popula- 
tions. As usual, the outlier is scenario E, which is bi- 
ased toward few cycles and long precession timescales, 



but all other scenarios produce 4.5 < 2V, 



cycles 



< 



and 



30 ms < (Tp rec (ti/2)) < 65 ms, for our fiducial a = 0.03 
case. We plot distributions of (N cyc \ cs ) in Fig. El where 
we see that low a values produce many more disk pre- 
cession cycles. 

In Figure [5] we plot contours to indicate the relative 
probability of BH-NS mergers producing combinations 
of T prec and ipd- These two quantities in general appear 
fairly uncorrelated, although in some scenarios (A, B and 
G) we do see a weak positive correlation, indicating that 
the most dramatically precessing SGRB disks will typical 
precess with longer periods. In no scenario do we see typ- 
ical T prcc (£i/2) values above 0.2 s; if we set aside scenario 
E, only a small fraction of events have T prec (i 1 / 2 ) above 
0.1 s. Scenarios C and H (soft NS equation of state) 
generally have the shortest precession periods, as a com- 
bination of rapidly spinning BHs and fairly low ip are 
required to disrupt these more compact NSs. Generally, 
?prec(ii/2) 0.01 s in all scenarios. 



V. DISCUSSION 

In this paper, we have explored much of the relevant 
parameter space for BH-NS mergers, and found that 
quantities relevant for disk precession are generally insen- 
sitive to the assumptions we have made about progenitor 
binary parameters. The important results of our cal- 
culations {(T pTec (t 1/2 )), (2V cyclos ), (ip d ), /grb) typically 
change by factors < 2 as we have varied our assumptions 
about the underlying populations of pre-merger BH-NS 
binaries. The one exception to this concerns pre-merger 
BH spin: in order for a significant fraction of BH-NS 
mergers to produce SGRBs, cibh <; 0.5 is required. Be- 
cause of the difficulty of increasing BH spin through pre- 
merger mass transfer |49j . this is equivalent to requiring 
modestly large natal obh values. The strongest preces- 
sion effects (i.e. large misalignment ipd and short period 
Tprec) arise from populations with top-heavy BH spin dis- 
tributions and stiff NS equations of state. Our results are 
notably insensitive to choice of BH mass function, so long 
as the mass function peaks near 6M Q , as is suggested by 
observations of X-ray binaries. 

On the other hand, our results do depend fairly sensi- 
tively on the details of how the disk will viscously spread 
outward. Using the a viscosity parametrization, we have 
seen that effective a values > 0.1 will strongly suppress 
the number of observable precession cycles, to the point 
where jet precession will rarely if ever be detectable. 
Likewise, a > 0.1 could be large enough to induce a 
Bardeen-Petterson warp in the thinner, inner regions of 
the disk, limiting global precession. We note again, how- 
ever, that what matters most for the time evolution of 
7p r cc is not the value of a in the inner regions of the 
disk, where GR effects enhance the turbulent viscosity 
produced by the MRI [36], but rather the outer regions 
of the disk, where local shearing box simulations that find 
lower (0.01 < a < 0.1) levels of MRI-generated viscosity 
are more appropriate [35) . 




FIG. 4. Probability distributions of N cyc \ es for scenarios A-H. As in Fig. 1, the dotted black lines represent a = 0.1, the dashed 
magenta lines a = 0.03, and the solid cyan lines a = 0.01. There is relatively little variation between progenitor scenarios, with 
the notable exception of scenario E (slow spins). The total number of precession cycles accumulated depends strongly on a. 



If jets align with the disk angular momentum axis, 
then they will precess around the total angular momen- 
tum vector by an angle w V'di because Jbh is signifi- 
cantly larger than Jdisk- In this case, observations of 
SGRBs associated with BH-NSs will often be marked 
by a clear "lighthouse effect," so long as Vd > 6>j Ct . 
This seems plausible, as observations of jet breaks in 
SGRBs suggest typical opening angles of ~ 10° [8J. If 
■0d < 0jct, then jet precession would, typically, be en- 
coded more subtly as a variation in the portion of the 
jet presented to the observer. Alternatively, if jets align 
with the BH spin axis, they will precess by a much 
smaller amount, since the angle between the BH spin 
vector and the total system angular momentum vector 
is cosi/'bh = (1 + j cos ipa) j^fl + 2jcosV>d + j 2 , where 
j = Jd/^BH- This angle will initially be smaller than tp^ 
by a factor ~ 10 if Jd ~ O.IJbh [22] , and this difference 
will increase slowly in time as angular momentum is ac- 
creted from the disk onto the BH (or quickly, if angular 
momentum is lost in outflows). Unless the typical 6?j e t is 
quite small, less than a few degrees, the dramatic light- 
housing effect previously discussed would be unlikely. A 
quasi-periodic signal will still be present in all these sce- 
narios, but it will be strongest for tightly collimated jets 
aligned with Jdisk- When we correct our results for the 



w 2 overestimate in ip^, we have found (ipd) ~ $jct, al- 
though there are substantial observational uncertainties 
regarding the distribution of jet opening angles. 

Detection of a precession signal appears observation- 
ally feasible: our study found a floor of ~ 25 ms for the 
"halfway" precession period (T prec (t 1 / 2 )}, but precession 
periods of ~ 50 msec were much more common, and in 
all cases the precession period grows in time. Both Swift 
and Fermi GBM have much finer intrinsic timing reso- 
lution, so the observability of precession signals will be 
ultimately limited by photon-counting statistics. Analy- 
sis of Fermi SGRB light curves finds evidence of structure 
in the ~ 10 — 50 msec range while Swift BAT also 
produces signals with structure in the 10s of msec |67| . 
For comparison, the BATSE archival data sample con- 
tains sufficiently high-resolution timing data to enable 
a search for ~ kHz frequency QPOs [SHIES]; indicating 
that the much longer precession signals should be clearly 
resolvable. 

The two primary challenges for observability are (1) 
the relatively low number of precession cycles, and (2) 
the rapidly growing precession period T prcc . However, we 
have shown that for most progenitor populations there 
exists a large-value tail to the N cyc \ es distribution, from 
which events with N cyc \ es w 5 can be observed, provided 



9 




FIG. 5. Each contour plot in this figure shows the range of "halfway" precession periods (vertical axis) and post-merger 
misalignment angles tpd (horizontal axis) for every BH-NS merger that produces a disk. T prcc is plotted in seconds, and in all 
scenarios is generally < 0.4s, although we note again that it will grow with time as the disk viscously spreads outward. As 
in previous plots, we show scenarios A-H. The logarithmic contour scale is included in panel (E); the range in colors covers a 
probability range of two orders of magnitude. In this plot we use the fiducial value a = 0.03. 



a < 0.1. Furthermore, searches for evolving quasiperi- 
odicity are feasible, provided the time axis of time series 
data can be rescaled to match appropriate theoretical 
models. In this paper we have employed a simple ana- 
lytic model for the viscous spreading of the disk, and for 
the first time have found that precession periods grow 
as i 4 / 3 . We hope that this will provide a starting point 
for searches for jet precession in SGRBs, but detailed 
hydrodynamical, and perhaps GRMHD, simulations are 
necessary to validate or refine this analytic expectation. 
We also note that both of these challenges are mitigated 
if a significant fraction of disk angular momentum is lost 
through an outflow, which will reduce the late-time scal- 
ing of T proc to i 4 / 5 . 

Pulsation of prompt emission may not be the only ob- 
servable implication of disk precession in SGRBs. A pre- 
cessing jet will sweep out a larger solid angle in the sky, 
enabling BH-NS mergers to make up a larger fraction 
of the observed SGRB rate than would be implied by a 
simple calculation (i.e. the intrinsic BH-NS merger rate 
multiplied by /grb)- This would likewise enhance ob- 
servability of BH-NS optical afterglows [70] . 

Unless there exists a large population of NS-NS bina- 



ries with millisecond spin periods (not accounted for by 
current population synthesis estimates), SGRBs due to 
BH-NS mergers will be distinguishable from those due 
to NS-NS mergers by the presence of a quasiperiodic sig- 
nal, with a typical period of ~ 30 — 100 ms. We have 
shown that this signal is robust to a large number of 
assumptions about the progenitor binaries, but evolves 
quickly and could become difficult to observe if jets align 
with the BH spin axis, or if the post-merger disk viscos- 
ity is large. A better understanding of both viscous disk 
spreading, and how GRB jet intensity varies with angle of 
observation, will aid future searches for this discriminant 
between SGRB progenitor binaries. 
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Appendix A: Innermost Stable Spherical Orbits 

We follow the simple formalism of Perez-Giz 2012 to 
calculate the ISSO radius for Kerr metric geodesies, and 
present it here for completeness (in this appendix we use 
geometrized units, G=c=l, for radial distance r and BH 
spin a). As mentioned in §2, the ISSO is the innermost 
stable orbit at constant radius, but fixed nonzero incli- 
nation around a Kerr BH. We define an inclination angle 
u such that C — cos l, the Carter constant Q = L 2 sin 2 t, 
and conserved vertical angular momentum L z = LC, i.e. 
l = corresponds to equatorial prograde orbits. Cal- 
culation of the equatorial plane ISCO, nsco, is well- 
documented in the literature [55] , and consists of finding 
the roots of the polynomial 

Z(r) = (r(r - 6)) 2 - a 2 (2r(3r + 14) - 9a 2 ) = 0, (Al) 

with one root the prograde and one the retrograde ISCO. 
The polar (t = ±7r/2) ISSO can be found at the root of 



P(r) =r A {r\r - 6) + a z (3r + 4)) 
+ a 4 (3r{r~2)+a 2 ) = 



(A2) 



that lies between r = 6 and r = 1 + \/3 + \/3 + 2^/3. 
Finally, the generic ISSO is a root of the polynomial 

S(r) = r 8 Z(r)+a 2 {l - C 2 )(a 2 (l - C 2 )Y(r) (A3) 
- 2r 4 A»), 

with auxiliary functions defined as 

X(r) —a 2 (a 2 (3a 2 + 4r(2r — 3)) (A4) 
+ r 2 (15r(r - 4) + 28)) - 6r 4 (r 2 - 4) 

Y(r) =a 4 (a 4 + r 2 (7r(3r - 4) + 36)) + 6r(r - 2) (A5) 
x (a 6 + 2r 3 (a 2 (3r + 2) + 3r 2 (r - 2))). 

Specifically, nsso is the root of S(r) located between the 
appropriate nsco (prograde or retrograde) and the polar 
ISSO. 

Appendix B: Post-Newtonian Merger Treatment 

The Ref. [5U] formulae discussed in §3 are presented 
in this appendix, along with a discussion of their ap- 
plicability to BH-NS mergers. These formulae are cal- 
ibrated to define the outcome of a BH-BH merger, al- 
though we apply them more generally to the case of 



BH-NS mergers with masses mi and rri2, and mass ra- 
tio q — mi/7712 < 1- We define the symmetric mass 
ratio ij = q/(l + q) 2 , and denote the dimensionless CO 

spin vectors <ij as having components aj and a 1 ", which 
are parallel to and perpendicular to the binary angular 
momentum, respectively. We further define 0± as the 
angle made between the radial direction and the vector 
A± = (mi + m2)(m2a,2 ± miai). Then the post-merger 
remnant mass is found to decrease by a fraction 



5M/M = T}E mC o + E 2 n 2 + E 3 n 3 



if 



(1 + 9) 2 

x (E s {4 + q 2 a ] l) + E A (1 - q)(a 2 l - qa 1 }) 

+E A \a 2 + gail 2 + E B \a£ + qa^\ 2 
x (cos 2 (0+ - e a ) + E c ) + E B \a 2 - gai| 2 

+E E \a 2 L - ga^| 2 (cos 2 (6_ - 3 ) + E F ) 



(Bl) 



Here the energy lost during the inspiral from infinity 
down to the plunge is fit as 



£1 



ISCO 



1 



0.103803?7 



1 



36^3(1 + q) 2 
5 



(B2) 



X (q(l + 2q)a)} + (2 + q)a l j) - 
\ y ' 1 1 ' 2 J 324^2(1 + q) 2 

x (a 2 . - 3(4 ) 2 - 2q{a 1 ■ a 2 - 344) 
+ 2 (a 2 -3( ai ) 2 )). 
The final spin vector of the post-merger BH is 



r-(, 



M J 



(B3) 



(Va(4 + 9 2 «i) + Jb(1 - <z)(4 - qa)l))n\\ 



+ (TTg) 2 

+ (1 - q)\a2 - ga ] L |v / ^ACOs(2(OA - 64)) + Jm"_l 
+ |4 + q 2 ai\VJs cos(2(O s - 6 5 )) + Ju S n± 



where the angular momentum radiated during the inspi- 
ral from infinity to the plunge is 



^isco 



2V3 - 1.525586277- 



1 



9^2(1 + q) 2 



(B4) 



(g(7 + 8q)a 1 } + (8 + lq)a. 



9^3(1 + a) 2 
x (aj - 3(4 f - 2q(a 1 ■ a 2 - 344) 



+^(^-3(4) 2 ) 



9^2(1 + q) 2 



{q{l+^ + {A + q)a 2 ) + l - a j^. 

The other variables in these formulae are empirical fit- 
ting constants, with values found to be E% = .341, E3 — 
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Source 


q 


xp 




"BH (NR) 


a' BH (PN) 


M d.NR 

M NS 


M d,F 
M-NS 


[55J 


1/7 


0° 


0.5 


0.67 


0.658 


< 0.4% 


0% 


55 


1/7 


0° 


0.7 


0.80 


0.786 


6% 


7.2% 


55. 


1/7 


0° 


0.9 


0.92 


0.913 


28% 


22.9% 


55 


1/5 


0° 


0.5 


0.71 


0.681 


6% 


6.5% 


32. 


1/3 


0° 


0.0 


0.56 


0.54 


5.2% 


5.21% 


32. 


1/3 


0° 


0.5 


0.77 


0.70 


15.5% 


16.3% 


32 


1 /Q 

1/3 





0.9 


0.93 


0.829 


38.9% 


28.3% 


32~ 


1/3 


20° 


0.5 


0.76 


0.699 


14.5% 


15.8% 


32. 


1/3 


40° 


0.5 


0.74 


0.691 


11.5% 


14.1% 


32 


1/3 


60° 


0.5 


0.71 


0.671 


8.0% 


11.4% 


.32. 


1/3 


80° 


0.5 


0.66 


0.636 


6.1% 


8.0% 



TABLE II. In this table we compare the full NR results for 
post-merger B H sp in a' B}i to the predictions of our PN fitting 
formula, Eq. (B3|. Here if) is the initial spin-orbit misalign- 



ment angle, q is the mass ratio, and the source column refers 
to the paper whose NR data we are comparing the PN results 
to. The agreement is quite strong for small q. For q = 1/3, 
roughly the upper limit considered in our BH mass functions, 
the error is < 10%, which is still acceptable for our purposes. 
In the last two columns we also demonstrate the good agree- 
ment between NR estimates for initial post-merger disk mass, 
Md,NR, and the prescription of Eq. |6|, Ma,f- 



0.522, E s = 0.673, E A = -0.3689, E A = -0.0136, £ B = 
0.045, £ c = 0,£ D = 0.2611, E E = 0.0959, £ F = 0, J 2 = 
-2.81, J 3 = 1.69, J A = -2.9667, J B = -1.7296, J A = 
Jm = Js = Jms = 0. We also set 62 = 63 = 64 = 65 = 
because of the weak dependence of the results on these 
parameters. With extremely low probability (~ 10 -5 ), 
Eq. (B3) can give super-extremal spin values, of a' > 1; 
we discard these cases from our sample when they ap- 
pear. 

As mentioned in §3, these formulae, which were de- 
rived, and calibrated, for BH-BH mergers, are found to 
give surprisingly good agreement with detailed results for 
BH-NS mergers simulated in full NR. We demonstrate 
the agreement in Table II, where we also plot the gener- 
ally good agreement between Eq. (|6| and NR results for 
the initial mass of the post-merger accretion disk. The 
largest errors in the disk mass fitting formula seem to 
occur for large disk masses (Md > 0.2Mns), and for high 
pre-merger spin-orbit misalignments (ip > 60° ) . The first 
of these cases occurs only for a small subset of our BH- 
NS mergers; the second occurs for a larger fraction, and 
may result in a modest overestimate of /grb- 



